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The quantum behavior of classically chaotic systems has been extensively studied both with time dependent and time 
independent Hamiltonians [TJ{7] . The main issue is that of determining fingerprints of classical chaos in the quantum 
Q_i mechanical behavior. For example, the spectral statistics of closed classically integrable [8 10_ and classically chaotic 
i 1 „ [TTl - fTB] quantum systems have been predicted to have clearly distinct properties. Many of the systems that are of 
physical interest are mixed, where some parts of the phase space are chaotic and some parts are regular. Spectral 
properties of mixed systems with time independent Hamiltonians were studied by Berry and Robnik |14| . In the 
present paper we study the classical/quantum correspondence properties of a mixed, open, time dependent system. 
(Here by "open" we mean that both position and momentum are unbounded.) 

The system we study consists of a particle kicked by a Gaussian potential defined by the Hamiltonian, 
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Models of this form were studied by Jensen who used it to investigate quantum effects on scattering in classically 
chaotic 1 15] and mixed |16| systems. This system can be experimentally approximated by a Gaussian laser beam 
acting on a cloud of cold atoms, somewhat similar to the realization of the kicked rotor by Raizen and coworkers |17| . 
As we will show, the study of Hamiltonian ([lj is particularly suited to the investigation of generic behavior of kicked, 
open, mixed-phase-space systems. In particular, we will focus on issues of fidelity [T5H2T| . decoherence [221 |2"3"] and 
scattering [2 M27|. Our main motivation in studying Hamiltonian (JTJ is that, with likely future technological advances 
(see Sec. | VII| for discussion), the phenomena we consider may soon become accessible to experimental investigation. 

Quantum mechanically, it is expected that classical phase space details on the scale of Planck's constant are washed 
out |281l2l)] . In contrast, one of our results will be that quantum dynamics can be sensitive to extremely fine structures 
in phase space, and this sensitivity is stable in the presence of noise [22, 23 . Phase space tunneling has been studied 
extensively [30-33 . For systems with many phase space structures complications arise due to transport between these 
structures. For our Hamiltonian (lj the motion is unbounded (i.e., the system is "open"), and therefore this system is 
ideal for the exploration of tunneling out of phase space structures and, in particular, for study of resonance assisted 
tunneling, a current active field of research [31-33 . 

The outline of our paper is as follows. Section |TT] presents and discusses our model system. Section [Hl| considers the 
quasi-energies of quantum states localized to island chains. Section [TV| introduces the fidelity concept and applies it 



to study different regions of the phase space including the main, central KAM island (Sec. IV A I, island chains (Sec 



IV B I, and chaotic regions (Sec. IV C I. Experimentally there is always some noise present in such systems. Also, 
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noise can be intentionally introduced. Section |V| considers this issue. Sectio n |VI| presents a study of the scattering 
properties of the system. Conclusions and further discussion are given in Sec. |VII| 

II. THE MODEL 
A. Formulation 

A particle kicked by a Gaussian beam is modeled by the Hamiltonian, Eq. (JTJ, with the classical equations of 
motion, 



p = -— = -K'T^e-^ Ht-Tn), (2) 



dx 

dH p 
dp m 

We rewrite the equations of motion in dimensionless form by defining x — x/A, i = t/T. The dimensionlcss momentum 
is correspondingly defined as p = pT / (mA). Thus we obtain the dimensionless equations of motion, 

oo 

p = -Kxe'^r S(i-n), (3) 

n— — oo 

X = p, 

where 

K'T 2 
m/X z 

Since in what follows we deal with the rescaled position, momentum and time we will drop the bar notation for 
convenience. By integrating ^ and defining p n = p{t = n._), x n = x (t = n._), where t = n_ is a time just before 
the n-th kick, we can rewrite the differential equations of the motion as a mapping, M, 

y[ . J Pn+i = Pn ~ Kx n e 2 y 
[Xn+l =Xn+Pn+l- 

The corresponding quantum dynamics in rescaled units is given by the Hamiltonian, 

2 00 

H = P --Ke-^ J2 Ht-n), (6) 
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where p — —ird Xl and r — HT/ (mA 2 ) is the rescaled h, namely, [x,p] — ir. The quantum evolution is given by 

irdtiff = -^d^if) - Ke~£ Ht-n)^, (7) 

or by the one kick propagator 

2 / K 

U\ = e _l fc exp f i — ) . (8) 

B. Properties of the Classical Map and the phase portrait 

In this subsection the classical properties of the map Eq. ^ will be presented. The first property is reflection 
symmetry, (x,p) — > (—x, —p). Phase portraits such as these presented in Fig. [ljand Fig. [2] are clearly seen to satisfy 
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Figure 1: (Color online) The phase space for K = 1. Colors (shades) distinguish different orbits. 




Figure 2: (Color online) The phase space for K — 4.5. Colors (shades) distinguish different orbits. 



this property. Like the standard map, Eq. ([5| can be written as a product of two involutions, M = J2</i, where 

Ji : (p,x) -> (-p,a; + p) (9) 
J2 : (p, x) — > p — Kxe~^~ , x^j . 

We will use this property for the calculation of the periodic orbits. From ^ we see that the only fixed point is 
(x = 0,p — 0). Linearizing around this point, we find that the trace of the tangent map is 2 — K. Therefore, this 
point is elliptic for < K < 4, and, for K = 1, the phase portrait of Fig. [T]is found, while for K > 4 this point is 
hyperbolic, leading to phase portraits like that of Fig. [2] Since the kicking as a function of x is bounded by K, for 
large initial momentum the particle is nearly not affected by the kicks, and continues to move in its initial direction. 
For < K < 4 we find a large island around the elliptic point (x,p) = (0, 0), and, for nearly all initial conditions near 
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x = p — 0, the motion is regular (i.e., lies on KAM surfaces). Further away from this fixed point, one finds island 
chains embedded in a chaotic strip. And even further away, the motion is unbounded. 



III. QUASI-ENERGIES OF AN ISLAND CHAIN 

In the semiclassical regime quasi-energies are related to classical structures. In this section we assume the existence 
of quasi-energy eigenfunctions u n (x). 

U lUn (x) = e- lE -u n (x), (10) 

such that u n (x) is strongly localized to an island chain of order r, and we attempt to calculate the quasi-energy E n . 
For this purpose we use the one-kick propagator U± to generate successive jumps in the island chain, 

Uiipi = ip i+ i, (11) 

where ipi is a wavefunction which is localized in island number i within the island chain. Further, we assume that 
this wavefunction can be expanded using the quasi-eigenstates of the island chain, 

tpi = / ] CjnUn jx) ■ (12) 
n 

Using this expansion we obtain a system of equations, 

U 1 ip i = ^c in u n (x)e~ lEn , (13) 



and 



Classically the i— th island is transformed to itself by r successive applications of the map M. In particular, an elliptic 
fixed point of the map M r is located in the center of the island. In the semiclassical limit the eigenstates of U[ are 
determined by M r and are close to the eigenstates of a harmonic oscillator centered on the fixed point of M r . The 
frequency of the oscillator, Vi, is such that the eigenvalues of the tangent map of M r , which transforms the i— th island 
to itself, are e Vi . This tangent map can be written in terms of the product of the tangent maps of M (i — > i + 1) , 
which transform the i— th island to the (i + 1) — th island. Consequently, since the eigenvalues are determined by 
the trace of the product of the tangent maps, they are independent of i (due to the invariance of the trace to cyclic 
permutations). In what follows we therefore drop the index i from vi. 
Choosing ipi as the eigenstate of U[, means that 

U[iPi = e^ip h (15) 
where j3 = tv/2 and we have taken ipi to be the ground state of the harmonic oscillator. Therefore, 

iPi = Y,Ci n u n {x)e-< E " r+ P). (16) 



Using the orthogonality of the u n (x), Eqs. (16 1 and (12 1 yield 



e 



The quasi-energies, obtained from (17 1 are, therefore, 



E n = — n + (3, 0<n<r, (18) 
r 

where /3 — —/3/r. Approximations to the quasi-energies can be calculated numerically by launching a wavepacket into 
one island in the island chain and propagating it in time, which gives 

iP(x,N) = Y,c ln u n (x)e- lE " N . (19) 



Taking a Fourier transform with respect to TV gives the quasi-energies. We have found that for K = 1 the chains with 
r = 8 and r = 16 accurately satisfy (18 1. 
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IV. FIDELITY 



The concept of quantum fidelity was introduced by Peres |18| as a fingerprint of classical chaos in quantum dynamics. 
It has subsequently been extensively utilized in theoretical |19| |20| [34J [35] and experimental studies [35H38] , for a 
review see [5T|. Most of this research has focused on the difference between chaotic and regular systems. Here we 
discuss fidelity for a mixed system. We have calculated the fidelity, 



S(t)= (Me lHlt/T e~ tH2t/T \M 



(20) 



where H\ 2 are Hamiltonians of the form (pi with with slightly different kicking strengths, ^1.2, and 4>o is the initial 
wavefunction. We note that the fidelity S (I) can be experimentally measured by the Ramsey method, as used in Ref. 
|37| . The fidelity is related to an integral over Wigner functions, 



S (t) = 2ttt J dxdpPfa (x,p) Pfc (x,p) , 



(21) 



where Pa, 



are the Wigner functions of cf>\2 = e lHl,2t ^ T <Po, respectively. 



We study separately the fidelity in the central island, in the island chain, and in the chaotic region (i.e., Eq.(20l 
with c/)q localized to these regions). 



A. Fidelity of a wavepacket in the central island 

First we prepare the initial wavefunction cf>o as a Gaussian wavepacket with a minimal uncertainty, namely, Ax 
Ap=(r/2) 1/2 , 



4>o (x) 



1 



2tt (Ax) 



1/1 



g-ipox/r exp 



(x - XqY 

4 (Ax) 2 



(22) 



We place cj>o in the center of the island, namely, xq = po = 0. Since the center of the wavepacket is initially at the 
fixed point, for Ax and Ap classically small, its dynamics are approximately determined by the tangent map of the 
fixed point. For this purpose we linearize the classical map ([5| around the point x = p = 0. This gives the equation 
for the deviations, 



The eigenvalues of this equation are, 



Sx n+1 



(1 - K) l\( Sx n 
-K 1 S Pn 



«1,2 



± 



(2 -JO ^iy/K(A-K) 



(23) 



(24) 



with 



u> = arctan 



VK (4 - K) 
(2-K) 



(25) 



which is the angular velocity of the points around the origin. In the vicinity of the fixed point, the system behaves 
like a harmonic oscillator with a frequency uj. Classically, the motion of the trajectories, starting near the elliptic 
fixed point, x = p = 0, stays there because the region is bounded by KAM curves that surround this point. For 
small effective Planck's constant, r, the quantum behavior is expected to mimic the classical behavior for a long time. 



Inspired by the relation between the fidelity and the Wigner function (see (21 1), we have defined a classical fidelity, 
S c (t), as the overlap between coarse-grained Liouville densities of Hi and H2 (this is similar to the classical fidelity 
defined in [39 ). To do this we first randomly generate a large number of initial classical positions using the initial 
distribution function, 



fo (x,p) 



1 

2?rAxAp 



exp 



' - -''(1 \" , fp-po^ ' 



A;i 



Ap 



(26) 
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Figure 3: (Color online) Quantum fidelity, S(t), (dashed red) and classical fidelity, S c (t) (solid blue). K\ — 1, Kq, 
t = 0.01, x = -0.25 and p = 0. 



1.01, 




Figure 4: (Color online) Classical density, which was initially placed at xo = —0.25 and po = after 500 kicks. Blue (dark) 
dots are for K\ — 1 and green (light) dots are for K2 = 1.01. 



corresponding to our initial 4>q given by ([22]). The coarse grained densities for Hi and Hi are then computed by first 
integrating these initial conditions and then coarse graining to a grid of squares in phase space of area r |40| . The 
motivation for this procedure is to check if structures in phase space of size smaller than t are of importance to the 
fidelity. A comparison between S (t) and S c (t) for xq — —0.25, po = 0, and t = 0.01 is presented in Fig. [3] The 
initial wavepacket is smeared on a ring in the phase space due to the twist property of the map. Since the probability 
density is preserved, the "whorl" which is formed contains very dense and thin tendrils. In Fig. [4] two such "whorls" 
are presented for Hi with K 1 = 1 and H2 with K2 = 1.01. When the two "whorls" coincide a fidelity revival is formed. 
Coarse graining the densities to a boxes of size r averages the differences between the two "whorls", obtained by Hi 
and Hi. This explains why the classical fidelity approaches 1 as the number of kicks becomes large. On the other 
hand, the quantum fidelity shows strong revivals which suggests that it feels the difference in trajectories between the 
two Hamiltonians. To understand the period of the revivals, we calculate, 8ui, the frequency difference between the 
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Figure 5: (Color online) A numerical (solid blue) and an analytical (dashed green) computation of the period of the fidelity 
revival as a function of K, 5K — 0.1, xq — po = 0, r = 0.01. 



two Hamiltonians, H\. 2 . Expanding uj around K\ gives 



u (K) = cj {Kx) + 



K — K\ 

y/K(A-K) 



0{{K-Kx) 



Therefore, the difference in angular velocity between two orbits of Hamiltonians, H\ and H2 is given by 

5ui = u{K 2 )-u(Kx) = 



^K 2 {A-K 2 ) 



(27) 



(28) 



for SK — K 2 — K\. This suggests that the fidelity, S (t) will be periodic, with the period T = tv/Suj. Note that we 
predict T = tt/Slj, rather than T = 2tt/Sui. This is because of the symmetry of the initial condition. Each point of Hi 
is chasing a point of H 2 which is its reflection through the origin of the phase space and, therefore, is found first at an 
angle of tt and not 2ir. To check this, we have calculated the period of the revivals numerically for < K < 4. First, 
fidelity was computed and Fourier transformed, then the second most significant value was taken as the period. In 
Fig. [5]we present a comparison of the analytic calculation of the period of the fidelity and the numerical computation. 
The correspondence is good through the whole range of the stochasticity parameter K but degrades near K = 4, 
where the elliptic point at the origin becomes unstable. Also, near K = 2, resonance chains appear near the fixed 
point x = p = 0, which results in poor agreement with the theoretical prediction, see Fig. [7] Very often it is assumed 
that quantum mechanical behavior is insensitive to phase space structures with areas smaller than Planck's constant, 
which results in an effective averaging on this scale |28l |2"5] . While this assumption is often correct [32], sometimes it 
is not |4lTl46] . The difference between S (t) and S c (t) demonstrated in Fig. [3] shows that fidelity may be sensitive to 
extremely small details in the classical phase space. In particular, a "whorP ^Sl affects the quantum dynamics. 
The small decay of the quantum fidelity seen in Fig. [3jis a result of tunneling. 

We stress that to observe the oscillations which appear on Fig. [3] requires sensitivity to the structure of the "whorl" 
of Fig. [4] In our quantum calculation the effective Planck's constant is r = 0.01 and it is obvious that the "whorl" of 
Fig. [4] exhibits structures on smaller scale, for example, in a square with sides of length 0.1 in phase space (of Fig. 
[4} one finds several stripes of the "whorl". Indeed, averaging over such a square leads to the classical fidelity that 
does not exhibit oscillations as the quantum fidelity does. We conclude that the structures on the scale smaller than 
the effective Planck's constant, t, are crucial for the oscillations in the quantum fidelity. Hence, structures of scales 
smaller than Planck's constant may dominate fidelity, which is a quantum quantity. 

For a wavepacket started around an initial point (xq,Po) 7^ (0,0) the behavior is similar, but with a slightly 
different period due to a decrease in the angular velocity for points far from the fixed point. Similarly to the case of 
(xq,Pq) = (0,0), we have calculated numerically the revival period for different values of K; this is shown in Fig. [6]. 
For K > 1.5 resonances appear near the launching point which introduce additional periods into the fidelity, making 
the analysis more complicated. 



8 



1400 



1200 - 



1000 - 



| 800 - 



600 - 



400 - 




200 1 1 1 

0.5 1 1.5 

K 



Figure 6: A numerical (blue circles) and an analytical (solid blue line) computation of the period of the fidelity revival as a 
function of K, 8K = 0.01, x = -0.25, p = 0, r = 2 x 10~ 4 . 
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Figure 7: (Color online) The phase space for K — 2.1. Colors (shades) distinguish different orbits. 

B. Fidelity for a wavepacket in an island chain 

We consider two different island chains occurring for different values of K. For K = 2.1 we have examined a 
chain of order r = 4 (see Fig. [7]) , and for K = 1 we have studied a chain of order r = 8 (see Fig. [T]). The initial 
wavepacket was launched inside one of the islands of the chain, and the we numerically computed the fidelity. In Figs. 
^{Kx = 2.10, K 2 = 2.11) and|9](i!fi = 1.00, K 2 = 1.01) we show the results of these computations. 

It is notable that there are three timescales in the graph of the fidelity. The shortest timescale is visible only in 
the inset of Fig. [9] and may be understood taking into account the symmetry of the equations of motion, x — > — x, 
p — > —p. This symmetry implies that each island has a "twin" which is found by reflection through the origin, 
x = p = 0. Therefore, the overlap between the islands of Hi and _ff 2 is a periodic function with a period of r /2, 
where r is the number of islands in the chain. Consequently, for the island chains used to obtain Fig. [8] and [9j the 
fidelity has periods of 2 and 4, respectively, on its shortest timescale. The intermediate timescale is due to a rotation 
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Figure 8: Fidelity of packet started inside an island chain of order 4. K\ — 2.10, ifjj = 2.11, r = 2 x 10 4 , and the center of 
the packet is started at x — 0.3198, p — 0, in the center of one of the islands of the chain. 
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Figure 9: Fidelity of packet started inside an island chain of order 8. K\ = 1, K2 = 1.01, r = 2 x 10 4 , and the center of the 
packet is started at x = 1.1312, p = 0, in the center of one of the islands of the chain. The inset is a zoom on the graph. 

of the wavepacket around the elliptic points of the island where it is initially launched. The central point in the 
island is a fixed point of M r . In r iterations, points in the island rotate with an angular velocity u)\ and ui2 for Hi 
and H 2 , respectively. The angular velocities can be calculated numerically by linearization of the tangent map of 
M r around the fixed point of the map M r . We find the fixed point by reducing M to a product of involutions pj, 
which allows us to reduce the search for the fixed points to the line p = in the phase space since any point on 
this line is a fixed point of J\ |47[ I48| . For K\ — \ and K2 — 1.01, the angular velocities are found to be uj\ = 1.10 
and 0J2 — 1.147. For K\ — 2.10 and K 2 = 2.11, the angular velocities are found to be lji = 0.391 and w 2 = 0.429. 
Therefore, the time it takes for a packet to accomplish a full revolution around the fixed points of M r is 2irr/uj, where 
u> = \ (uj\ + L02) ~ Wi ~ L02 (see Table M. The longest timescale of the fidelity is the timescale when the difference 
between the angular velocities is resolved T = 2nr/Suj. In Table [I] we compare those periods deduced directly from 
Fig. [8] and Fig. [9] and the periods calculated by finding w li2 from the tangent map. We see that the agreement is 
excellent. 
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longest period 
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Table I: This table compares two ways of calculating the periods of revivals for the resonance chains. In one way we have 
deduced them from the Figures [8|9| and in the other way we have calculated them using the tangent map. This is done for two 
different resonances: r = 4, for K\ = 2.1, Ki = 2.11; and r — 8 for Ki = 1, Ki = 1.01. For both cases r = 2 x 10~ 4 . 
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Figure 10: Fidelity of packet started inside a chaotic layer. K\ — 1, Ki = 1.01, r = 2 x 10 4 and the center of the packet is 
started at x — —2, p — 0. 



C. Fidelity of the wavepacket in the chaotic strip 



For the fidelity of a packet started inside the chaotic strip (see Fig 
which is dependent on K\. This is half a period in this chain/strip, 
which is a characteristic of chaotic regions 



10 1, we notice a strong revival after 6 kicks 



After this revival the fidelity decays to zero, 
Detailed exploration of this region is left for further studies. 



V. DEPHASING 



We now investigate the effect of dephasing by adding temporal noise to the time between the kicks. The classical 
equations of motion with the dephasing are given by 

Pn+i = Pn - Kx n e ^ ; (29) 

Xn+l = X n + (1 + Stn) ■ Pn+l, 

and the quantum one kick propagator is 

= e -*£(i+«») exp e-*£\ , (30) 

where St n is a random variable which is normally distributed with zero mean and a standard deviation <j t - The 
standard deviation of the St n , corresponds to the strength of the noise. We find that the noise results in an escape 
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Figure 11: (Color online) Quantum fidelity, S (t), (dashed red) and classical fidelity, S c (t) (solid blue) for a dephasing noise of 
strength a t = 0.01, K x = 1, K 2 = 1.01, r = 0.01, x = -0.25 and Po = 0. 




Figure 12: (Color online) Quantum fidelity, S (t), (dashed red) and classical fidelity, S c (t) (solid blue) for a dephasing noise of 
strength a t = 0.001, K x = 1, K 2 = 1.01, r = 0.01, x = -0.25 and p = 0. 



outside of the island, which yields additional decay in the fidelity. Since we are interested in the difference between 
the two wavefunctions only inside the main island, for each kick we normalize the wavefunctions of Hi and Hi such 
that their norm is equal to 1 inside a region of \x\ < Xb — 3. This gives the following expression for the fidelity 



S(t) = 



-iHxt/r 



0o (x')) (e-Wrfa (x')) dx> 



1/2 ' 



with he classical fidelity S c (t) defined in a similar way. We have numerically calculated the fidelity for the same 
situation as in Fig. [3] with added relative noise of a t = 0.01 (Fig. 11) and a t — 0.001 (Fig. 12 1. We notice that the 
noise introduces additional decay in the quantum fidelity. 
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Figure 13: (Color online) Quantum fidelity, S(t), (dashed light red) and classical fidelity, S c (t) (solid dark blue) for two 
different realizations of a dephasing noise of strength a t = 0.001, K\ = K-z = 1, t = 0.01, x — —0.25 and p = 0. 
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Figure 14: (Color online) Classical density, which was initially placed at xq = —0.25 and po = after 5 x 10 4 kicks, K\ — K2 = 1. 
Colors (shades) correspond to two different realizations of a dephasing noise of strength at = 0.01. 



To isolate the effect of noise from the decay in the fidelity due to the difference between Ki and K 2 we set K\ = K 2 
and use two different noise realizations with the same strength a t ■ From Fig. [13] we notice that classical fidelity 
initially decays very fast due to the noise and than slowly recovers approaching a value of 0.8. This is due to the 



coarse graining to the scale of r. To illustrate this we plot in Fig. [14] the classical densities after 5 x 10 4 kicks for a 
packet initially launched at xq = —0.25. We notice that the densities for the two Hamiltonians highly overlap, which 
explains the high fidelity. In Fig. [15] we observe the corresponding quantum wavepackets. Contrary to the classical 
fidelity, the quantum fidelity decays rather slowly with the noise, suggesting that it is more robust to noise than the 
classical fidelity. 
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Figure 15: (Color online) Wavepackets, which were initially placed at xo = —0.25 and po — after 5 x 10 4 kicks, K\ — K2 = 1. 
Colors (shades) correspond to two different realizations of a dephasing noise of strength at = 0.01. 



VI. SCATTERING 



We now investigate the difference between quantum and classical scattering behavior by studying the evolution of 



a wavepacket initialized outside of the main island of the phase space, Eq. (22 1 with xq = —2, po = 0, Ax = Ap — 
(t/2) 1 / 2 . In the classical case both the classical chaos, as well as the numerous small island structures introduce, 
an erratic behavior for the transmission and reflection coefficients as a function of the initial launching position and 
energy [T51 Due to effective phase space smoothing of areas much smaller than our effective Planck's constant, 
r, we expect that fine scale fractal-like features in the classically erratic scattering dependence will be averaged out. 
To quantify this behavior, we measure the transmission and reflection coefficients for a wavepacket defined as the 
transfered or reflected probability mass, either quantum or classical. Classically, it is the fraction of initial trajectories 



(generated using ( 26 1 ) reflected or transmitted by the main island for a given time, while quantum mechanically, we 



measure the total escaped probability up to time t from the island area, |x| < Xb, 

ft I— Xb 



I) 



L(t) = I dt' I \ip (x,f)f dx (31) 



R(t) = / dt' / |V>0M')| 2 dx, 

JO Jx b 

where Xb is the margin of the main island (we choose Xb — 4) , L (t) and R (t) are probabilities to be scattered to 
the left or the right of the island till time t, correspondingly. To determine those probabilities, we use the continuity 
equation for the probability, 

d t \ ( \^\ 2 dx\=Tlm[^d x iP*)\ x =b-^d x r)\x=a], (32) 



so that, 



L{t) = jj dt 1 I dt" {^d x r -rd x i>)\ x= ^ Xb , (33) 



Jo 

t 



R(t) = ~ I dt'(ipd x ip*-rd x Tp)\ x = 







2i 

In Figs. |16||21| we compare the quantum and classical scattering of a wavepacket launched from the left of the 
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Figure 16: (Color online) Total quantum (solid blue line) and classical (blue dots) probabilities for scattering to the left of the 
island [x < —xi) as a function of the number of kicks. K = 1, t = 0.01, xq — —2, po = 0. 
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Figure 17: (Color online) Total quantum (solid blue line) and classical (blue dots) probabilities for scattering to the right of 
the island (x > Xb) as a function of the number of kicks. K = 1, r = 0.01, xq = —2, po = 0. 



main island. We notice that there is a substantial difference, which decreases when we decrease the effective Planck's 
constant, r. Figures 16]|T8 and Figs. T&pT differ in the initial launching position of the wavepacket (xq = —2. for 



Figs. T&fl8 and xq = —3 for Figs. I&pT l. We notice that the scattering is sensitive to xq. Different aspects of chaotic 
scattering for this problem were explored in [16] , and in particular, the effect of small h on washing out rainbow 
singularities of the classical scattering function. 
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Figure 18: Total quantum (solid blue line) and classical (blue dots) probabilities to stay in the island (\x\ < Xb) as a function 
of the number of kicks. K = 1, r = 0.01, x = —2, p — 0. 
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Figure 19: Total quantum (solid blue line) and classical (blue dots) probabilities for scattering to the left of the island (x < —Xi) 
as a function of the number of kicks. K = 1, r = 0.01, xq = —3, po = 0. 



VII. DISCUSSION AND CONCLUSIONS 

A. Discussion of experimental realizability 

In the present work the classical and quantum dynamics of a system with a mixed phase space were studied. It is 
proposed to realize this system by injecting cold atoms into a coherent, pulsed, gaussian light beam. The phase space 
structures, which can be seen on Figs. Tj2 and [7] are controlled by the parameters of the beam via the parameter K. 



Since it is relatively straightforward to control the parameters of gaussian beams, the proposed system is ideal for 
the exploration of dynamics of mixed systems. In what follows limitations on experimental realizations are discussed. 
First we consider the realizability of an approximately one dimensional situation necessary for the validity of our 
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Figure 20: Total quantum (solid blue line) and classical (blue dots) probabilities for scattering to the right of the island (x > Xb) 
as a function of the number of kicks. K = 1, r = 0.01, xo = —3, po = 0. 
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Figure 21: Total quantum (solid blue line) and classical (blue dots) probabilities to stay in the island 
of the number of kicks. K = 1, t = 0.01, xq = —3, po = 0. 



la; I < Xb) as a function 



theoretical results. Let us assume that the gaussian beam propagates in the z direction. Its profile in the xy plane is 



2A2 2A?, 



(34) 



Assuming that the extent of the light beam is much smaller than the Rayleigh length, zr = ttA 2 /X where A is 
the wavelength, and the z dependence of the potential can be ignored. The potential of Eq. (34 1 can be well 
approximated by exp (— x 2 / (2A 2 )) in ([lj, for sufficiently small values of y 2 /A 2 , and, to facilitate this, it is appropriate 
to consider A y 3> A, i.e., a quasi-sheet-like beam. Such beams arc experimentally realizable via routine methods. 
To analyze this situation, the normalized map M of (|5| should be replaced by one with exp (— x 2 /2) replaced by 
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cxp x 2 /2 + y 2 (A/Aj,) 2 j. In addition, there are equations for y n and its conjugate momentum p Vt n, which in 
dimcnsionless units with y and p y<n rescaled by A and T / (mA), respectively, take the form, 



Vn+1 — Vn + Py, n+1, 



2 ~2 ST V-n 



(35) 



where 



Ky=K 



A 
AT 



(36) 



Since K « 1 and A/A y <C 1, it can be assumed that K y <C 1. Therefore, the motion in the y direction is slow relative 
to the motion in the x direction. Thus exp (— £ 2 /2) can be approximated by its time average (exp (— x 2 / (2A 2 ))) = p 
(which is of order unity), and, for sufficiently small y the y— motion (351 can be described by a Harmonic oscillator 
with a force constant 2K y p <C 1. Conservation of energy E y implies that the maximal value of y satisfies 



E,. 



(37) 

The energy E y is determined by the initial preparation. Let us assume that initially the atoms form a Bose-Einstein 
Condensate (BEC) and are in a harmonic trap that is anisotropic where the frequency in the y direction is v[ y in 
experimental units, and v y — Tv' y in our rescaled units. We assume that the center of this trap y satisfies y ^ Umax- 
The experiment starts when the trap is turned off. Assuming the atoms are in the ground state, their energy in our 
rescales units is 

T 2 



\ hv 'y 



iA 2 



= ^VyT < Ey. 



(38) 



We desire the effect of the motion in the y direction on the motion in the x direction (Eq.Q with exp (— x 2 / (2A 2 )) 
replaced by V of ( 34 1) to be negligible. Thus it is required that 

v = yl^(^) «1- (39) 



In this case the y motion corresponds to a variation in K of the order AK • 
reduces to 



Kr\. Using (37) and (38 1, condition (39 1 



\VyT_ 

4 K 



< 



E v 
2K 



(40) 



where, since we are interested only in crude estimates, we have replaced p by one. The initial spread in y is given by 
the ground state of the harmonic oscillator, where (y 2 ) = t j {2v y ), and we require that the expectation value of y 2 
satisfies (y 2 ) <§C j/^ ax , resulting in 



A 
A^ 



2v, 



< T}. 



For both inequalities ( 40 ) and (41 1 to be satisfied it is required that 

2 



2,/ 



< v. 



< V- 

T 



The resulting fundamental lower bound on i] is 



T 

8K 



< rj. 



(41) 



(42) 



(43) 



Reasonable experimental values are A/A y « 10 2 and v y k, 0.1. For r = 10 2 and K « 1 the lower bound on rj 



For v. 



is 10 5 leaving a wide range for 'engineering' of BEC traps so that the v y is in the range ( |42 1 
t = 10~ 2 and K « 1 we can make r\ < 10~ 3 . Since this value of rj is small compared to the value of AK 



s 0.1 and 

used in our fidelity calculations (Figs [3] [5] [TT] [l2| ) , those calculations are expected to be uneffected by y motion for 
our assumed parameters. It is also encouraging to see that noise of a higher level does not destroy fidelity oscillations 
(see Fig. 12 1. One should note, however, that the variation of the effective K of the motion in the x direction is 
slow, with effective frequency A/ A y that for A/ A y « 10~ 2 is of order 10~ 2 . For these reasons, we expect that, the 
model that we have explored theoretically in the present work should be realizable for a wide range of experimental 
parameters. 



18 



B. Conclusions 

The main result of this paper is that the quantum fidelity is sensitive to the phase space details that are finer than 
Planck's constant, contrary to expectations of Refs. [28, 29 . In particular, the fidelity was studied and predicted to 
oscillate with frequencies that can be predicted from classical considerations. This behavior is characteristic of regular 
regions. Fidelity exhibits a periodic sequence of peaks. For wavepackets in the main island, it was checked that the peak 
structure is stable in the presence of external noise but the amplitude decays with time. For wavepackets initialized 
in a chain of regular islands, it was found that the fidelity exhibits several time scales that can be predicted from 
classical considerations. For wavepackets initialized in the chaotic region, the fidelity is found to decay exponentially 
as expected. It was shown how quasi-energies are related to classical structures in phase space. Substantial deviation 
between quantum and classical scattering was found. These quantum mechanical effects can be measured with kicked 
gaussian beams as demonstrated in the present work. 
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